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1 Analytic Perturbation Theory in QCD 



Abstract 

We give a short introduction to the Analytic Perturbation Theory (APT) [TJ [51 [31 HI [S] 
- - - , in QCD, describe its problems and suggest as a tool for their resolution the Fractional 

^ ■ APT (FAPT) [11|71[H1IS]. We also describe shortly how to treat heavy-quark thresholds in 

I FAPT. As an applications of this technique we discuss (i) the pion form factor calculation 

' in the Euclidean FAPT and (ii) the Higgs boson decay bb in Minkowskian FAPT. We 

conclude with comparison of both approaches, Euclidean and Minkowskian FAPT. 

O 

O 

(N . 

In the standard QCD Perturbation Theory (PT) we know that the Renormalization Group (RG) 
' (-h'I equation das[L]/dL = —a^ — . ■ ■ for the effective coupling UsiQ"^) = a[L]/Pf with L = ln((5^/A^), 
CL|' f3f = bo{Nf)/{47r) = (11 — 2N j / 3) / {Airj^ Then the one-loop solution generates Landau pole 
p ^. singularity, a(i)[L] = 1/L with subscript (/) meaning /-loop order. 

D I In the Analytic Perturbation Theory (APT) we have different effective couplings in Minkow- 

skian (Radyushkin [10] and Krasnikov&Pivovarov [11]) and Euclidean (Shirkov&Solovtsov [3]) 
regions. In the Euclidean domain, —q^ = Q^, L = InQ^/A^, APT generates one set of images for 
the effective coupling and its n-th powers, whereas in the Minkowskian domain, 

1/^ I = s, Lg = Ins/A^, it generates another set, {2l„[Ls]}„gpj. APT is based on the RG and 

causality that guaranties standard perturbative UV asymptotics and spectral properties. Power 
^ , series '^^'^(i) [-^1 transforms into non-power series ^m^^m^^il^] APT. 

By the analytization in APT for an observable /(Q^) we mean the "Kallen-Lehman" rep- 
resentation 

^ . with the spectral density defined through the perturbative result, Pf{cr) = (l/7r)Im[/P*'''*(— a)] . 
^ I This results in different analytic images in Euclidean and Minkowski regions 

An[L] ^ AeK[L]] = r^lM^da, (2a) 

Jo cr + Q'^ 

2l„[L,] ^ AmK[L]] = r^da. (2b) 

Js 0- 

Then in the one- loop approximation spectral density is p^\a) = 1/ [ln^(o"/A^) + vr^] and ana- 
lytic couplings are 



*e-mail:bakulev@theor. jinr.ru 

^We use notations f{Q^) and f[L] in order to specify what arguments we mean — squared momentum or 
its logarithm L = ln{Q^ / A^), that is f[L] — /(A^ ■ e^) and is usually referred to Nf — 3 region. 
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whereas spectral density for the analytic images of the higher powers (n > 2,n e N) is p^n\s 
(l/(n — 1)!) {—d/dLs)^~^ P^i\s) and analytic images of the higher powers are: 



1 



(n- 1)! 



d 

dL 



^ Ai[L] 



(4) 



2 Fractional APT 

In the standard QCD PT we have also: 

• the factorization procedure in QCD that gives rise to the appearance of logarithmic factors 
of the type: a''[L]L-M 

• the RG evolution that generates evolution factors of the type: B{Q^) = \^Z {Q"^) / Z 
B{ii'^), which reduce in the one-loop approximation to Z{Q'^) ~ a'^[L] with = 7o/(26o) 
being a fractional number. 

That means that in order to analytize perturbative QCD expressions we need to construct 
analytic images of new functions: a'^ , o'^L™',.... The result of this procedure is just the 
Fractional APT (FAPT). Formally speaking, in order to construct FAPT we need to use formulas 
([2]) with substituting spectral densities with integer indices n, namely Pn(c), by corresponding 
spectral densities Pv{(y) with fractional indices v. But in the one-loop approximation FAPT 
can more effectively be constructed using recursive relations (jll), see for details in [6, 8j. The 
explicit expressions for couplings in Euclidean, ^[/^^[L], and Minkowskian, 2l[/^^[L], domains are 
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F{e-^,l- v) 



sm 



{u — 1) arccos iL/^Jn'^ + L"^ 



(5) 



Here F{z, v) is reduced Lerch transcendental function, which is an analytic function in u. 
Interesting to note that ^[.^^[L] appears to be an entire function in v, whereas 2li^'*[L] is de- 
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Figure 1: Graphics of (left panel) and 21^^'' [L] (right panel) for fractional v e [2,3]. 

termined completely in terms of elementary functions. These couplings have the following 
properties: 



(i)r 



(6a) 



^« [L] ^ for ^ . N ; 2lL^i[L] ^ L , ^ - ^ , ^U^-.^L,...^ (6b) 
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A^^^l 



-L] 
-L] 



for m > 2 , m S N ; 



(6c) 



^First indication that a special "analytization" procedure is needed to handle these logarithmic terms appeared 
in [12] , where it has been suggested that one should demand the analyticity of the partonic amplitude as a whole. 
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[iboo] = 21^^) [±oo] = 0form>2, mGN; 
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We display graphics of A^u"^ [L] and 21^^^ [L] in Fig. [D one can see here a kind of distorting mirror 
on both panels. Next, in Fig. [2]we show graphics for = 2, 3, 4, 5. Here we can trace the partial 
values 



0; 



(7a) 
(7b) 



Graphics for w4i/[L] as functions of v at fixed values of L can be found in our last papers [B]. 
We compare the basic ingredients of FAPT in Table [T] with their counterparts in conventional 
PT and APT. 
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Inverse powers (al^])^" — ^_m[L] = L™ 2l_„[L] = + 0(^2) 



Index derivative a'^ln^a — P^'A 



Logarithms a'^L'^ — ^i/-fc 2l^-fc 



Table 1: Comparison of PT, APT, and FAPT in Euclidean (FAPT(E), L = In (QV^^)), and Minkowski 
(FAPT(M), L ^ hi (s/A^)) domains in the one-loop approximation. In the row, named 'Inverse powers', 
we put = L'" + 0(7r2) that encodes just Eq. (|gbl) . 



3 Developments of FAPT 
3.1 Two-loop coupling 

In this section we want to show how good is FAPT in approximating the two-loop analytic cou- 
pling expressions by expanding them in non-power series in terms of one-loop analytic couplings. 
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To this end, we remind the two-loop equation for the normahzed couphng a(2) = 6oa/(47r): 



da{2) [L] 
dL 



2 

-0(2) 



[L] [l + ci 0(2) [L]] with ci = -i . 



(8) 



The exact solution a(2) [L] of this equation satisfies the following functional equation: 



0(2) m 



+ ci In 



a(2) [L] 



1 + cia(2)[L] 



L. 



Its exact solution is known [i3\ to be 



cil + W.i{zw[L]) ' 



(9) 



(10) 



where zw[L] = (1/ci) exp (— 1 + ivr — L/ci) and Wk, k = 0, ±1,..., denote different branches 
of the Lambert function W{z), defined through functional equation z = W{z) e^^^\ Review of 
its properties can be found in [TH [T31 US] • 

To construct FAPT for the two-loop quantities one should use formulas ([2]) with substituting 
the spectral densities Pn{<y) by the corresponding two-loop spectral densities pi {a), defined by 



-Im 



TT 



a^2) [L - ivr] 



fill 



Then analytic images in Euclidean, [L], and Minkowskian, domains are defined 

through 

We can also expand 0(2) [i] in terms of 0(i)[L] = \/L with inclusion of terms 0{a?^-^-^): 

0(2) [L] = 0(i) [L] + ci oji) [L] In a(i) [L] + c\ \L\ (in^ 0(i) \L\ + In 0(i) [L] - l) + . . . (13) 
and then produce analytic version of this expansion 

4^)^"^^^[L] = 4^) + ciP422 + c?(l^^ + P-l)4i + ---- (14) 



(2)r 



da . 



a 



(12) 
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Figure 3: Left panel: Comparison of relative errors A2"^^"'"[L] (dotted line) and Af'^^"'"[L] (solid line) in 
FAPT. Right panel: Comparison of relative errors AP[L] (dashed line) in standard PT and Af^P'T[L] 
(solid line) in FAPT. 
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In Fig. [3] we demonstrate nice convergence of this expansion using relative errors of the 2- and 
3-term approximations: 

.4S"[Z,]+ciPyl!.'ijL] 



aFAPT 
^2 



[L] 



1 



A 



FAPT 



L] 



A FAPT 
^2 



1 



a(i) [L] + ci a2^) [L] In a(i) [L] + cf a^^^ [L] (in^ a(i) [L] + In a(i) [L] - 1 

a(2) W 



(15) 
(16) 
(17) 



We see that relative accuracy of the 3-term approximation in FAPT (see the left panel of Fig. 
E]) is better than 2% for L > —2. In the same time, the right panel of Fig. [3] demonstrates 
that relative accuracy of the same 3-term approximation in standard PT even at L ~ 1 is much 
higher — about 10%, whereas in FAPT it is smaller than 1%! 
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Figure 4: Left panel: The solid line corresponds to yl2^^[L], computed analytically via Eq. (fT8|) : dashed 
line represents the result of a numerical integration, while stars correspond to the available numerical 
results of Magradze in [15]. Right panel: The solid line represents -4o.62 [^]' computed analytically 
via Eq. (fTH|) . while the dashed line stands for -'^■62 r 



We obtain the corresponding expansion for the two-loop coupling with index v. 



FAPT 



^P^+P-l 



A^^Ul] + .... (18) 



and display comparison of different results for ^2 [-^] on the left panel of Fig. HI On the 

right panel of this figure we show comparison of FAPT and standard QCD PT with respect to 
the fractional index (power) of the coupling, fixed at the value v = 0.62, corresponding to the 
evolution exponent of the second moment of the pion distribution amplitude. 

In the Minkowskian region convergence of the FAPT expansion for the two-loop coupling 



2^(2);FAPT ^ ^(1) ^2ciV Sl^g [L] + c? [3 p2 + 2 p - 2] [L] + 



(19) 



is also very good, but in the vicinity of the point L = (Landau pole in the standard PT) it 
is not so fast, so that we need to take into account 0(c^)-terms in order to reach 5% level of 
accuracy, for more details look in [6]. 

Now we turn to the corresponding expansion for the two-loop coupling with logarithm: 



<Jm=AE[(a(,))^L] =4'ii + ciP4^)-c?4'|i + f 4?2-^-4S3 + 0(c?) . (20) 



1(2) 



(2) 



l(2) 



,(2) 



The exact spectral density can be easily found, 

^(2) 



-^(2) [^] 



(211 
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with R(^i^2)[L] and (/?(i 2)[-^] being inverse modula and phases of the corresponding one- and 
two-loop densities. In Fig. [5] we show the relative deviations 

, , 41i + 0(ci)+Q(cf) + 0(cf) + 0(cf) 

^3,41-^1. 31, ij = i • [^^) 

-^1.31,1 

Once again, we see that relative accuracy of the 0{cf) approximation in FAPT is better than 
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Figure 5: The solid line corresponds to ^2^''[i], computed analytically via Eq. (fT5|) : dashed line repre- 
sents the result of a numerical integration. 

2% for all L with the largest deviations in the vicinity of the point L = 0. In the Minkowski 

(2) 

region convergence of the FAPT expansion for the corresponding quantity is also very 

good, but not so fast as in the Euclidean case, so that we need to take into account 0(cf )-term 
in order to reach the same level of accuracy, for more details see in [D]. 




3.2 Heavy-quark thresholds 

Construction of FAPT with fixed number of quark flavors, Nf, is a two-step procedure: we start 
with the perturbative result [a(Q^)] , generate the spectral density Pi/{(t) = (l/7r)Im [a'^(— o")] ., 
and then obtain analytic couplings Au[L] and via Eqs. ©• Here A'^^ is fixed and A'^^- 

dependent coefficient 6o(A^/) is factorized out. We can proceed in the same manner for Nf- 
dependent quantities: [ag{Q^; Nf)Y =^ Pui'^'^^f) = Pu[Lcj]Nf] = Pu{o')/Pf ^ -^A^',]^/] and 
2lj^[L; Nf] — here Nf is fixed, but not factorized out. 

Global version of FAPT, which takes into account heavy-quark thresholds, is constructed 
along the same lines but starting from global perturbative coupling [a|'°''(Q^)] being a con- 
tinuous function of due to choosing different values of QCD scales Aj, corresponding to 
different values of Nf. We illustrate here the case of only one heavy-quark threshold at s = ml, 
corresponding to the transition Nf = 3 ^ Nf = 4. Then we obtain the discontinuous spectral 
density 

p^°^ia) = pt^'iLa] =e{K< U) Pn [W, 3] + 9 (L4 < La) Pn [La + A4; 4] , (23) 

with La = In ((t/A|), Lf = In ^my/A|^ and A/ = In ^Ag/A^^ for / = 4, which is expressed in 

terms of fixed-flavor spectral densities with 3 and 4 flavors, Pn[L] 3] and pn[L + A4; 4]. However 
it generates the continuous Minkowskian coupling 

^t'^'lLs] = ^(L,<L4)(2t,[L,;3]-2l,[L4;3]+2t,[L4 + A4;4]) 

+ e{U<Ls) 2l,[L, + A4;4]. (24) 

and the analytic Euclidean coupling (for more detail see in [9]) 

U _ 

Arm = ML + M-A]+ / ^"'-[^-^j-p- [L +A4;4] ^^^^ ^^^^ 
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We analyze now how important is the deviation of the global FAPT from fixed Nf FAPT. 
In the Euclidean domain we may write 



^f^lL] =A^[L + A4; 4] + AA^[L] 



with 



L4 



' 1 + e^'^"^ 



(26) 



(27) 



whereas in the Minkowski domain — 

Slf^iL] = 214L+A4; 4] + A2l^[L] (28) 

with 

A2l^[L] = 0(L < L4) [214L;3] -2l^[L4;3] +2l^[L4 + A4;4] -2l^[L+A4;4]l (29) 
In Fig. [6] we show the relative values of deviations in Euclidean and Minkowski domains. We 
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Figure 6: Relative deviations of global coupling with respect to fixed- A^/ coupling in APT. Solid lines 
correspond to the ratios A^^i [i] /^f ""^ [L] (leftt panel) and A2li [L]/2lf (right panel). 

see that in both domains deviations vary from —20% at large values of —L ~ 10, in the vicinity 
of -L ~ —5 they go through the zero, then grow up to the value +20% which is reached at L ~ 0, 
and then tends to at L ^ 00. 



4 Electromagnetic pion form factor at NLO 

The scaled hard-scattering amplitude truncated at the next-to-leading order (NLO) and evalu- 
ated at renormalization scale /i^ = XrQ^ reads [THl [IZl [HI [12] 



as {XrQ^) 



47rQ2 



4')(x,y;^2^/Q2) 



(30) 



with shorthand notation {x = 1 — x 

fWf^ .,.,,2 /^2n _^_,(0) 



t\^> {x, y; fi'p/Q') = Cf C i^^ v) \^ + In {xy)) ^ + ^0 4''^^ ^r) + *H i^^v) ■ (31) 
The leading twist-2 pion distribution amplitude (DA) [20] at normalization scale is given 

by m 



,(FG), 



ip^{x, fi'j,) =6x{l-x) 1 + 02(^1.) C^^"i2x - 1) + aiifij,) C2'"i2x - 1) + . . . 



2 X r^3/2. 



(32) 
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All nonperturbative information is encapsulated in Gegenbauer coefficients anifJ^'p)- 

To obtain factorized part of pion form factor (FF) one needs to convolute the pion DA with 
the hard-scattering amplitude: 

FJ-*(Q2) = ^^(^; ^ T^LO y. q2^ ^ ^^(y. . (33) 

X y 

In order to obtain the analytic expression for the pion FF at NLO in [22] the so-called "Naive 

(2) 

Analytization" has been suggested. It uses analytic image only for coupling itself, A\ , but 
not for its powers. In contrast and in full accord with the APT ideology the receipt of "Maxi- 
mal Analytization" has been proposed recently in [23]. The corresponding expressions for the 
analytized hard amplitudes read as follows: 



In Fig. [7] we show the predictions for the factorized pion FF in the standard pQCD and 




[GeV'] [GeV'] [GeV'] 

Figure 7: Left panel: Factorized pion FF in the standard MS scheme. The solid line corresponds to 
the scale setting 11^ — 1 GeV^, dashed line — to fi^ = Q^, dotted lines — to the principle of minimal 
sensitivity, whereas dash-dotted lines — to the principle of fastest apparent convergency. Central panel: 
Factorized pion FF in the "Naive Analytization" . Right panel: Factorized pion FF in the "Maximal 
Analytization" . On both panels solid lines correspond to the scale setting /ip = 1 GeV^ , dashed lines — 
to /i^ — Q"^, dotted lines — to the BLM [21] prescription, whereas dash-dotted lines — to the ay-scheme. 



in "Naive" and in the "Maximal Analytization" approaches. We see that in the "Maximal 
Analytization" approach the obtained results are practically insensitive to the renormalization 
scheme and scale-setting choice (already at the NLO level). 

We show also the graphics for the whole pion FF, obtained in APT with the "Maximally 
Analytic" procedure using the Ward identity to match the non-factorized and factorized parts 
of the pion FF, see the right panel of Fig. [7| The green strip in this figure contains both 
nonperturbative uncertainties from nonlocal QCD sum rules [25\ [77] and renormalization 
scheme and scale ambiguities at the level of the NLO accuracy. 

It is interesting to note here that the FAPT approach, used in [7j for analytization of the 
ln((5^//Li|.)-terms in the hard amplitude ()30p . diminishes also the dependence on the factorization 
scale setting in the interval //|. = 1 — 50 GeV^. 

This conclusion starts to be even more pronounced in the complete FAPT analysis of the 
factorized pion FF. Indeed, if we put fi^ = then we obtain in the pion FF convolutions with 
(pTr{x,Q'^) which contains ERBL evolution factors a2n(^o) [q^s(Q^)/q^s(a*o)]'^^" with z^2n(-^/) = 
7o(2n)/(26o(A^/). Numerically U2 = 0.62-0.72 and 1/4 = 0.90-1.06. In T^LO [x,y,Q'^) we have 
three types of contributions: as(Q^), a1{Q'^), and bQ{Nf)a^{Q'^). The scheme of "analytization" 
of A'^y-dependent quantity is the same as in the case of the global approach, see in Fig. [8l Indeed, 
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we start with bo{Nf)ag{Q'^), obtain corresponding spectral density pf.";,^!^), being in case of the 
single heavy-quark threshold of the type 



' (L, < L4) 60(3) Pn [La, 3] +9 (L4 < La) bo(4) Pn [La + A4; 4] 



(36) 



and then construct corresponding analytic image A^°^^ [L] , using integral representations ([2]) 
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Figure 8: Schemes of "analytization" of A^y-dependent quantities. On the left panel the starting object 
to analytize is [af°^{Q'^)]'', whereas on the right panel — bQ{Nf)[af°^{Q'^)]'' . 

In Fig. [9] we show the predictions for the factorized pion FF in the complete FAPT ap- 
proaches with different settings for //p (see figure caption for details) and for comparison — 
in the "Maximal Analytization" scheme with pp = Q^. We see that both approaches produce 
practically the same results. 
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Figure 9: Factorized pion FF in the complete FAPT approach. The solid line corresponds to the 
factorization scale /i| ~ Q^, dashed line — to = 6 GeV^, whereas dotted line — to the "Maximal 
Analytization" with /i| = Q^. 

We conclude this section with the note that in the Euclidean domain taking thresholds 
and A^j-dependence of coefficients into account generates tiny correction! The main advantage 
of the complete FAPT approach: There are no problems with thresholds. Bonus: Pion FF 
automatically appears to be analytic function out of the Minkowski cut f9]. Note here that 
in this case one should use in the Minkowski region not the coupling 2ti^(s), but instead the 
analytic image of the Euclidean coupling on the lower side of the cut, namely, Ai,{—s — iO) = 
ILeAu{s) + i-Kpy{s). In contrast to the real coupling 21,^(5), this quantity naturally generates 
an imaginary part of the pion form factor, which is measured in e^e~ tt~^tt~ reactions. 
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5 Higgs boson decay into bb pair in FAPT 

In this section we analyze the Higgs boson decay to a bb pair. For the decay width we have 



(Mh) 



Gp 
4^/27r 



MnRs{Ml) with Rs{Ml)=ml{Ml)Rs{Ml) 



(37) 



and Rs{s) is the i?-ratio for the scalar correlator, see for details in (HJ [28]. The running mass 
m{Q'^) is described in the three-loop approximation by the renormgroup equation [9j 



[as{Q^)Y' [1 - ^21 asY'' [1 + 522 ol^' , (38) 
8 GeV) and i/q = 1-04, ^21 = 0.672, 1^21 = 



with RG- invariant mass ttl^^-^ (for 6-quark mf,.(3) 
-0.743, 622 = 0.029, and 1^22 = 8.59 at Nf = 5. 

In the standard PT direct multi-loop calculations are usually performed in the Euclidean 
region for the corresponding Adler function Da, = 3 [l + X^„>o ^^n a" (Q^)] j where QCD per- 
turbation theory works. Functions Ds and Rs{s) = 3 [l -|- X]„>o^" ctsl'S)] can be related to 
each other via a dispersion relation without any reference to perturbation theory. This gen- 
erates relations between r„ and dn coefficients, involving the famous vr-^-terms due to integral 
transformation of ln'^((5^//i^) in dn- 



The influence of these vr^-terms can be substantial, as has been shown by Baikovei al. |28j . 



5.1 Comparing different approaches to calculate rfj__>fefe(AfH) 
We compare now the results of different approaches to calculate i?g(M^). 



Baikov et al. (BCK) [28] used standard QCD PT at the 0(af)-order with Aj^Lg 
231 MeV: ^ 



R. 



(i=4)BCK 



3m 



(i=4) 



a. 



(i=4) ■ 



n>l 



vr 



(39) 



Broadhurst et al. (BKM) [29j in the approach of the so-called "Naive Non-Abelianization" 
(NNA) used the "contour-improved" optimization of power expansion. Their formulas are 
very closed to the one- loop FAPT (see more detailed discussion in [U [30] ) : 



0-2 
^^{1=1) 



1 ^ai;,glob, s 

n>l 



(40) 



(1)'Z 

We show the prediction obtained in this approach with using A)^'Z^ = 111 MeV (extracted 

from condition 21^^^''''°'' (m|) = 0.120) by the dotted ( green) line in Fig. 1101 It appears to 
be lower than the standard PT result (dashed red line) by ~ 8%. 



In the FAPT approach with "freezed" at the value dn{5) coefficients dn{Nf) (as also has 
been done in the BCK and BKM approaches) we define ^n+uoi^) ^ ttis analytic images 
of al'~^''° {Q'^)[l + 61 asiQ'^)Y^ t^e Minkowski region and obtain 



m(3);gioi 



3 n>(3);giob/ 

■(.) + E^^(5)^^ 

n>l 



(41) 
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Figure 10: Results of calculations of rjj^5f,(My) in different approaches. On both panels dashed (red) 
lines display the BCK results and dotted (green) lines — the BKM results. Solid line on the left panel 
corresponds to the results obtained in the FAPT approach with freezed coefficients c?„ = dn{Nf = 5), 
whereas on the right panel — to the results of the complete FAPT approach. 

(3) 

The corresponding predictions for the value Ajy^^^ = 220 MeV (being in accord with the 
last analysis of precision electroweak data by the LEP Collaboration [31J and normalized 
to 2l(^^^''^'°'^(m|) = 0.118) are represented in the left panel of Fig. [IO]by the solid (blue) 
line. We see that it is very close to the BCK curve (difference is of the order of 1.5%, and 
as compared with three-loop BCK result — even of the order of 1%). 

• In the complete FAPT approach with the complete analytization of A^j-dependencies we 
obtain 



' n>l 



TT" 



(42) 



where the analytic images *Bn+i/o;d„(s) absorb all A'j-dependence of dn coefficients. Here 
^7v]=5 ~ '^'^'^ MeV in order to match R{m\) value by the LEP Collaboration [31] in 
the complete FAPT approach: for this reason 

2^{3);giob^^|^ = 0.119 and differs slightly 

from the "freezed" FAPT value 21^^'*''^'°'" (m|) = 0.118. The results of this approach are 
shown in the right panel of Fig. [IHlby the solid (blue) line and appear to be smaller than 
the "freezed" results. This difference varies from 14% (at Mh = 50 GeV) to 18% (at 
Mh = 150 GeV). 

We can conclude from this comparison that the standard pQCD power series and the FAPT 
non-power series expansions appears to be quite close to each other in the region L ^ 1 
in the scenario with "freezed" coefficients dn{Nf), corresponding to the value Nf = 5. The 
complete FAPT setup, which bears in mind the whole Nf dependence of the perturbative results, 
naturally generates another, in our case smaller, result. This is not a surprise, because such an 
analytization approach effectively average coefficients dn{Nf) in the regions with different Nf 
values, resulting in reductions of their values and, hence, of the whole non-power series sum. In 
order to make this effect more transparent we define effective coefficients, corresponding to the 
complete FAPT series i?g^^'^^^'^[L] as follows: 

Sg{3);giob ry^i 

df[L] = "+;°^""^ \ (43) 

In Table [21 we show these effective coefficients: one can immediately realize that they indeed di- 
minished by 16 to 18%. Physically this effect of the complete analytization of Rs{s) corresponds 
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Table 2: The effective coefficients d^[L], see ()43p . at L = 11 — 13. In two last columns we show 
the corresponding values of Rg [L] in cases of the "freezed" and the complete FAPT. 







n = 


77/ 1 


n = 2 


71 

/ V O 


p(3);FAPT(5)r.n 
rig [jjj 


^(3);FAPT. . 
rig [Jvj 


dn{Nf = 


= 5) 


1.00 


5.67 


42.0 


353 






df[L = 


10] 


0.85 


4.93 


36.4 


296 


32.44 


27.78 


df[L = 


11] 


0.83 


4.80 


34.9 


276 


29.11 


24.41 


df[L = 


12] 


0.82 


4.66 


33.2 


252 


26.38 


21.64 



to taking into account contributions of loops with t-quarks even in the region -^s < 175 GeV, 
where in the standard pQCD in the MS-scheme only quarks u, d, s, c, and b contribute. We 
see that this effect is sizable. 

6 Concluding Remarks 

We conclude with the following resume: 

• The implementation of the analyticity concept (the dispersion relations) from the level of 
the coupling and its powers to the level of QCD amplitudes as a whole generates extension 
of the APT to FAPT. 

• The rules how to apply FAPT at the two- and three-loop levels are formulated. 

• The convergence of the perturbative expansion is significantly improved when using non- 
power FAPT expansion; 

• We formulate the rules how to account for heavy-quark thresholds in FAPT. 

• As an additional advantage we obtain the minimal sensitivity to both the renormalization 
and factorization scale setting, revealed on the example of the pion electromagnetic form 
factor. Threshold problem (how to fix the value of 6o(-^/) in the 0(Q;g)-term) resolved. 
The result of the complete FAPT prescription appears to be very close to the result of 
the Brodsky-Lepage-Mackenzie prescription, see in Figs. [7] and M 

• In application to decay H'^ bb we revealed that the complete FAPT analytization 
prescription reduces the results by ~ 16% as compared with the "freezed" FAPT ones. 

The comparison of FAPT in Euclidean and Minkowski regions, done in Sections 2 and 3, 
shows us close analogy between the two regions — the effect, named by Shirkov and Solovtsov as 
effect of "distorting mirror" . But the results of the complete FAPT applications to the analysis 
of the pion form factor in the Euclidean domain (Section H]) and of the Higgs boson decay width 
in the Minkowski domain (Sectional) are formally quite different: For the pion form factor case 
the complete FAPT prescription generates approximately the same result as standard FAPT 
with "freezed" A^j-dependent terms [23^, 21, see in Fig. El whereas for the Higgs boson decay 
this approach generates reduction of the order of 16%. 

In order to explain these difference we analyzed the Adler function Ds{Q'^) with the same 
perturbative coefficients dn{Nf). Then we revealed that the magnitude of the corresponding 
reduction of the complete FAPT results in the Euclidean region in comparison with the "freezed" 
FAPT ones is of the same order (18%) as in the Minkowski case. This demonstration proves 
that the magnitude of this reduction is related with the strong dependence of perturbative 
coefficients dn{Nf) on Nf and not with specific character of analytization procedure in the 
Minkowski region. 
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